Ridge Formation Induced by Jets in pp Collisions at 7 TeV 
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An interpretation of the ridge phenomenon found in pp collisions at 7 TeV is given in terms of 
enhancement of soft partons due to energy loss of semihard jets. A description of ridge formation 
in nuclear collisions can directly be extended to pp collisions, since hydrodynamics is not used, 
and azimuthal anisotropy is generated by semihard scattering. The observed ridge structure is 
then understood as a manifestation of soft-soft transverse correlation induced by semihard partons 
without long-range longitudinal correlation. Both the pr and multiplicity dependencies are well 
reproduced. Some predictions are made about other observables. 

PACS numbers: 25.75.Gz, 13.85.Ni 



I. INTRODUCTION 

The observation of ridge structure in two-particle cor- 
relation in pp collisions at 7 TeV by the CMS Collabo- 
ration at Large Hadron Collider (LHC) [l| has opened 
up the question of whether it has a similar origin as that 
already found at Relativistic Heavy-Ion Collider (RHIC) 
in Au-Au collisions at 0.2 TeV 0tQ|- A great deal is 
known about the ridge in heavy-ion collisions, since vari- 
ous experiments have studied two-particle (with or with- 
out trigger) and three-particle correlations. The domi- 
nant theme is that the ridge exhibits the effect of high 
or intermediate-pT jets on a dense medium. If the phe- 
nomenon seen at LHC reveals similar features upon fur- 
ther investigation, it would imply that soft partons of 
high density can be created in pp collisions and can af- 
fect the passage of hard partons through them. If not, a 
new mechanism needs to be found. Various theoretical 
speculations have been advanced with varying degrees of 
attention to the specifics of the CMS data [S-Q- I n this 
article we propose a model that is an extension of our 
past interpretation of the ridge phenomena in the RHIC 
data, but is particularly suitable for pp collisions at LHC, 
since the dynamical origin is jet production rather than 
hydrodynamics. We have a simple formula that can re- 
produce the CMS data quantitatively with the use of two 
parameters that can clearly describe the physics involved. 

The most direct approach to the study of ridges is 
to consider only events selected by triggers with p^ 18 in 
an intermediate px range, as first reported by Putschke 
3, 10]. The dependence of the ridge yield on centrality in 
nuclear collisions indicates that the ridge is formed when 
there is a jet in a dense medium. Having an exponen- 
tial behavior in p^ soc at values less than p^ Ig suggests 
that the ridge particles are related to the soft partons, 
but they have an inverse slope larger than that of the 
inclusive distribution, implying an enhancement effect of 
the jet [1, [H|. If triggers are not used as in the study of 
autocorrelation, ridges are also observed at IA77I > 1 in 
central collisions @,|B]. For pp collisions at LHC we can- 



not presume the existence of a dense medium of partons, 
which is a possibility we leave open. However, we can 
and shall assume that ridge formation is due to high- or 
intermediate-£»T jets, whether or not the jets are detected 
by triggers. Our goal is to study the properties of cor- 
relation generated by semihard jets. It should be noted 
that there are models in which the ridge phenomenon can 
occur without jets, such as in Refs. |9l. [12h15| . 



In the hadronization model based on Refs. [TT|, the 
ridge component (due to the recombination of thermal 
partons) manifests the effect of the semihard parton on 
the medium. The soft partons have exponential depen- 
dence on the transverse momentum kr, whose inverse 
slope is T in the absence of semihard partons. For the 
ridge component the inverse slope is increased to T" > T 
due to the enhancement of the thermal motion of the soft 
partons caused by the energy loss of the semihard par- 
ton that passes through the medium in the vicinity [l Tt j . 
That is soft-semihard correlation, which we shall apply 
to even pp collisions where the notion of thermal partons 
may be questionable. It is known empirically that there 
exists an exponential peak at small pr at LHC [18l420l ]: 
that is sufficient for us to refer to the underlying par- 
tons as soft, the recombination of which gives the \ow-px 
hadrons. 



In Sec. II we give a short summary of our past work on 
ridges with emphasis on the distinction between trans- 
verse and longitudinal correlations. It is significant to 
note that the data on ridge reported by PHOBOS Q do 
not imply the existence of long-range longitudinal corre- 
lation upon closer examination. In Sec. Ill the transverse 
correlation is extended to |A?7| > 1 appropriate for CMS 
measurement. Quantitative analysis of the ridge yield in 
pp collisions is then carried out in Sec. IV. In the last sec- 
tion we give the conclusion along with some predictions. 
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II. TRANSVERSE AND LONGITUDINAL 
CORRELATIONS 



Longitudinal correlation has been the prim ary concern 
of most theoretical studies on ridges [UH, HJH^ . The 
observation by PHOBOS [1] that |A?y| can be as large 
as 4 has led to the conclusion that there is empirical 
evidence for long-range correlation, which is an inherent 
property of flux-tube models. There are, however, two 
other aspects about the ridge structure that one should 
also consider in addition to the large- A77 aspect of the 
PHOBOS data. One is A: the property of ridge in the 
small Aj] limit, and the other is B: the question of how 
large should A77 be in order for the correlation to be 
regarded as long-range. We comment on them in the 
context of what have been observed at RHIC as a prelude 
to our discussion about the ridge found at LHC. 

A. Transverse Correlation. At midrapidity dihadron cor- 
relations in the azimuthal angles have been studied in 
detail at RHIC; in particular, the dependence of the az- 
imuthal correlation on the trigger angle <f) s relative to 
the reaction plane reveals features that are important 
about ridge formation [24T - I27T ] . Any model on the origin 
of ridges at \An\ > 1 should contain properties that are 
consistent with the azimuthal behavior at \Arj\ < 1, since 
all observed ridge structure have common behavior in Acj) 
throughout the An range. The ridge yield as a function of 
4> s has been studied in a model where the angular correla- 
tion between the trigger and local flow direction is limited 
[28l | . It is found that a Gaussian width of a ~ 0.34 can 
reproduce the data [24| [26|, H3] ■ The model suggests that 
thermal activities of the soft partons in the vicinity of 
the trajectory of the semihard parton (i.e., within a cone 
of angular range of a) are enhanced by the energy loss 
of the latter to the medium. Those enhanced thermal 
partons hadronize into the ridge particles that rise above 
the background. That is transverse correlation between 
the soft and semihard partons, the only type that can be 
studied when \Arj\ is restricted to < 1. After finding sat- 
isfactory explanation of the azimuthal correlation in the 
data this way for triggered events, the natural question 
to follow is how such correlation influences the single- 
particle distribution when triggers are not used. Semi- 
hard partons can be pervasive if their kx is around 3 
GeV/c or lower. It is found that the semihard-soft trans- 
verse correlation can give rise to a significant azimuthal 
anisotropy (UHIl; and that V2{pr, N valt ) can be quanti- 
tatively reproduced as a consequence of the ridge effect in 
inclusive distribution |30T |. This will become a key input 
in our discussion below where the nature of the transverse 
correlation will be made explicit. 

B. Longitudinal Correlation. At first sight of the PHO- 
BOS data on the A77 range of the ridge distribution 
anyone having some familiarity with multiparticle pro- 
duction is likely to regard 1 772 — 771 1 ~ 4 as indicative of 
long-range correlation between the trigger at r\\ and ridge 
particle at 772. However, to quantify the notion of corre- 
lation range it is important to compare it to the ?7-range 



of the single-particle distribution. A recent study shows 
that the ridge distribution in A77, denoted by dN^ 1 /dAn, 
can be related empirically to the single-particle distribu- 
tion, dN^Jdmby using the two relevant sets of PHOBOS 
data only [J, [31| without any theoretical input |32| . That 
phenomenological relationship 
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ch 



dm 



(1) 



involves a shift in 772 of the charge hadron and an inte- 
gration over the trigger 771 , and shows that the range of 
correlation in A77 is no more than the range of the inclu- 
sive distribution apart from the smearing of the trigger 
acceptance, which lengthens the A77 range by 1.5. The 
implication is that there is no long-range longitudinal cor- 
relation. Any successful model of ridge formation should 
be able to explain the simple relationship shown in Eq. 
(Q]) . In Ref. [32j an interpretation of that relationship is 
given in terms of transverse correlation that we discuss 
in more detail in the next section. 



III. RIDGE AT IA77I > 1 

The phenomenological verification of Eq. ([1]) directs 
one's attention to the origin of ridge formation without 
intrinsic longitudinal correlation at large A77. From all 
that have been learned experimentally about the ridges, 
there is no indication that such structure can be found 
in the absence of any jet. Even in autocorrelation stud- 
ies where no triggers are used, ridges are found in the 
kinematical region where minijets are detected 0]. Our 
approach is therefore to start with jet-induced transverse 
correlation at \An\ < 1 and to extend it to larger 77 sep- 
aration, in contrast to other studies where long-range 
longitudinal correlation at low px exists without jets 
and then a large-pr parton is introduced to define the 
A77 range. The approach we adopt was actually advo- 
cated even before the discovery of ridge was reported by 
Putschke [Io| at a time when the phenomenon was re- 
garded as the pedestal lying under the jet peak |33l. 1 341) . 
Now, with more data and model analyses of the trans- 
verse correlation at hand, the extension to large A 77 can 
be done with more definiteness. 

To be more specific, let us consider the ridge found by 
CMS at LHC, where only charged particles with \rj\ < 2.4 
and pt < 4 GeV/c are used in the analysis. In that ac- 
ceptance region the hadron p^ is less than 22 GeV/c, so 
Feynman xf is < 6.3 x 10~ 3 at ^/s — 7 TeV, and the 
corresponding partons that recombine have even lower 
x values. Those are soft wee partons deep in the sea, 
whose correlations can be strongly influenced by fluctua- 
tions. Suppose that a semihard scattering occurs in a pp 
collision at 7 TeV and sends a parton to the 77 ~ region 
with a parton momentum kx in the 5-10 GeV/c range, 
which we shall regard as intermediate at LHC. Whatever 
the medium effect on it may be, it can lead to a cluster 
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of hadrons with limited range in r\ and <j) [lj . It cannot 
directly cause the production of an associated particle at 
rj = 2.4 since the pl of that particle can exceed 20 GeV/c, 
hence forbidden by energy conservation. Any particle 
produced outside the jet peak carries longitudinal mo- 
mentum that is driven by the initial partons (right- or 
left-movers) of the incident protons. In the conventional 
parton model it is assumed that there are no significant 
longitudinal constraints on those initial partons [H, ■ 
We add, however, that their transverse momentum distri- 
bution can be affected by the semihard scattering before 
they recede from one another. At early time the right- 
and left-movers need not be arranged as in Hubble ex- 
pansion, i.e., a right-moving parton may be located on 
the left side of the region of uncertainty, and vice-versa; 
hence, those initial partons can be sensitive the passage 
of the semi-hard parton across their ways. The quantum 
fluctuations that generate the transverse kx distribution 
of the forward (or backward) moving partons may be en- 
hanced by the energy loss of the semihard parton. More 
specifically, let exp(— Ht/T) represent the kx distribution 
in the absence of semihard scattering; then our asser- 
tion is that the distribution changes to exp(— kr/T') with 
T" > T in the presence of semihard scattering, provided 
that the affected partons are in the vicinity of the semi- 
hard parton trajectory in the transverse plane, i.e., A<j) 
is limited on the near side. Furthermore, such a change 
occurs for all partons independent of their longitudinal 
momenta up to x ~ 10~ 2 , say. This enhancement is in 
essence the transverse correlation discussed in Sec. II. A, 
but now the semihard parton at r\ sa induces a change 
in the transverse distribution of the soft partons from T 
to T' at all rj in the limited region \rf\ < 2.4 under study. 

The CMS experiment does not identify any particle as 
the trigger, so the pseudorapidity of the semihard par- 
ton cannot be specified. All charged particles accepted 
in the window \r)\ < 2.4 are used for the analysis of the 
two-particle correlation. Thus the correlated particles 
may be at 771.2 = ±2.4, resulting in |A^| = |r?i — 772 j 
as large as 4.8. Hereafter, r/i and r] 2 do not refer to 
trigger and associated particles, respectively, but to any 
two particles whose correlation is measured by CMS. The 
semihard parton may be anywhere in between ±2.4. The 
huge jet peak observed in Ref. [l[ corresponds to particles 
that are produced by thermal-shower recombination and 
therefore must be close in rj to the semihard parton, but 
the peak distribution in A77 does not indicate where it is. 
The flat ridge distribution that lies below the jet peak 
only reveals the response of the medium in terms of en- 
hanced thermal partons without any information about 
the locations of the shower partons. The ridge particles 
have transverse distribution that is characterized by the 
same inverse slope T" as for the enhanced soft partons. 
That is a property of recombination [l6|, [3(| • No explicit 
longitudinal correlation has been put in. 

In order to describe pion and proton production in the 
same formalism of recombination of thermal partons at 
low px, it is shown that the replacement of pr by Et, 



where E{pt) = {rn^+p^) 1 / 2 —rrih, h — it or p, is sufficient 
to account for the mass effect and that the inclusive ridge 
distribution can reproduce 1$ C®ri -W por t) at low Et [301 ) . 
Being the difference between the enhanced distribution 
and the background, that ridge distribution is 



i?(p T )=i?o(e-^ /T ' 



-E T /T\ 



(2) 



for nuclear collisions. It is the soft response to the semi- 
hard partons. We will apply the same description to pp 
collision below. The difference AT = T' — T is a measure 
of the magnitude of the influence by semihard scattering 
without which there is no ridge. 



IV. RIDGE YIELD IN pp COLLISION AT LHC 

We now focus on the ridge yield measured by 
CMS. Let the single-particle distribution be p(pr, v) — 
dN/pTdr/dpT, which will be abbreviated by pi(i) for 
the iih particle, so that two-particle distribution is de- 
noted by p2(l)2). Define two-particle correlation by 
C 2 (l,2) = p 2 (l,2) - |C?i(i)pi(2). The measure for ridge 
used by CMS is 



ftcMs(l,2) = ArC 2 (l,2)/pi(l)pi(2), 



(3) 



where N is the number of charged particles in a multi- 
plicity bin. In more detail the quantities in Eq. ([3]) are 
averaged over bins of pt, so Ref. [l[ exhibits 



ft C Ms(PT,A?7,A(/>) = N- 
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J [pt] dp Tl PTiPi{i) 





(4) 



where [pt] denotes the range of integration from px — 0.5 
to p T + 0.5 (GeV/c). A ridge then appears in the 2D A77- 
A<p distribution. A projection of it onto A</> is done by 
integrating | A77I over the range 2.0 to 4.8. The associated 
yield in the ridge is then determined by integrating over a 
range of Acj> around where 7?-cms is above its minimum, 
i.e., 



Y R (p T ,N) 



dA(f) 



±4.1 



dArj K C ms{pt, At7,A</>).(5) 
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This measure of the ridge yield is given for 4 bins of pt 
and N each The data points are shown in Fig. I. 

What is remarkable about the data is that Yr is very 
small for both 0.1 < pt < 1 and 3 < pt < 4 GeV/c, 
but jumps up by nearly an order of magnitude in the 
1 < pt < 2 GeV/c bin. It is very unusual in high-energy 
physics where the pt behavior is so drastically different 
on the two sides of 1 GeV/c. The increase of Yr with N 
is not surprising, especially if one has in mind that jets 
are connected with the ridge phenomenon. 

Our explanation of the pt and N dependencies of Yr 
is very simple, based on what has already been discussed. 
We assume no longitudinal correlation, as in [32], which 
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FIG. 1: Ridge yield vs multiplicity N for 4 bins of pr- Data 
are from Ref. and lines are from model calculation. 



can explain the PHOBOS data Q. Thus the only con- 
tribution to 6*2(1,2) is from transverse correlation that 
gives rise to the ridge distribution given in Eq. ^ as an 
^-independent response to the semihard jet at any r]j e t- 
Wc therefore write 



C 2 {1,2) = R(1)R(2). 



(6) 



This is a very unconventional description of correlation 
that we are proposing, since one usually expects an un- 
factorizable form for correlation. The two particles at r\\ 
and 772 are correlated because their pt distributions are 
both enhanced by the jet. R(l) and R(2) are indepen- 
dent responses, so they enter into C2(l,2) as factorized 
products. We emphasize that Eq. ^ is a correlation be- 
tween two soft particles, each of which being correlated 
transversely to the unobserved jet as described by Eq. 

An analogy for this is the adage that rising tide 
raises all boats — even though, we add, there are no in- 
trinsic horizontal correlations among the boats. Putting 
Eq. JBJ in (gj and (0 we obtain 



Y R (p T ,N) = cNf[ 



J [pt] dpriPTiR(PTi,N) 



/ brl dp T iPTiPl{pTi) 



(7) 



where c is an adjustable parameter that depends on the 
experiment. This is an explicit formula that enables us 
to do phenomenological analysis. 

The single-particle distribution for \rj\ < 2.4 at 7 TeV 
is given by CMS in the Tsallis parametrization 



Pi{pt) = Po(l + -=r) 
nl 



(8) 



with To = 0.145 GeV/c and n — 6.6. The average pr 
found from the above fit is (pr) — 0.545 GeV/c. 

We use Eq. ([5]) in (0 and fit the data in Fig. 1 with two 
parameters (apart from normalization) , which we choose 
to be T and /?, where 



AT 



PkxN, 



AT = T' — T. 



(9) 



This dependence on N is reasonable, since at higher N 
there is higher probability for jet production and hence 
larger AT, which is in the exponent in Eq. ([2]). The result 
of the fit is shown by the solid lines in Fig. 1 for 



T = 0.294 GeV and (3 = 0.0175. 



(10) 



Evidently, our model reproduces the data very well for all 
Pt and N bins. Yr{pt-, N) is small at small pr because 
R(pt) in Eq. © is suppressed as pr — > 0. The reason for 
that is discussed below. Yr(pt,N) is also small at large 
Pt', that is due both to the exponential suppression of 
R(pt) and the power-law decrease of piipr) at high pt- 
The increase with N that is most pronounced in the 1 < 
Pt < 2 GeV/c bin, where R(pt) is maximum, is clearly 
due to the enhancement of T when jet production is more 
likely in accordance to Eq. ©. At N = 100, AT/T is 
about 8%, which is slightly lower than that observed in 
nuclear collisions at RHIC where T — 355 ±6 MeV/c and 
V = 416 ± 22 MeV/c for 4 < p T rig < 6 GeV/c Q. 

The reason why R(pt) must vanish as pt — > is re- 
lated to azimuthal anisotropy in nuclear collisions. We 
have advocated the view that the ridge component before 
being averaged over <j) contains all the dependence of 
the inclusive distribution fl7l W\ . In that approach which 
has been worked out in more detail recently in [301 ] . it is 
shown without using hydrodynamics that v 2 (referred to 
as elliptic flow in hydro description) can be reproduced 
at all centralities, provided that R(pt) — > at vanishing 
Pt because v 2 (pt) — > 0. Since the azimuthal behavior is 
determined primarily by the initial geometry of the col- 
lision system [T3, 0, H3| , such an approach may well be 
applicable to pp collisions, for which the validity of hy- 
drodynamics used for nuclear collisions is doubtful. The 
origin of the <p dependence in the geometrical approach 
is the anisotropy of semihard emission when the initial 
configuration is almond-shaped. Similarly, it is reason- 
able to consider the initial configuration in pp collisions 
also, when the impact parameter is non-zero, and we ex- 
pect significant 4> anisotropy in the produced particles. 

The Tsallis distribution in Eq. ((SJ) has the property 
of a power-law behavior at large pt, but an exponential 
behavior, exp(— Et/To), at low pt- It is then of interest 
to note the difference between the values of T and T, 
the latter being twice larger than the former. It may 
appear as being inconsistent; however, the average (j>t) 
of exp(-£r/T) is 0.6 GeV/c, only 10% higher than that 
for Eq. ((SJ). Thus different parametrizations of the Et 
distribution give essentially the same physical quantity. 
Eq. dSJ) is a fit of the CMS data PU that emphasizes the 
p^jJ 1 behavior at high pt, while Eq. @ is a theoretical 
model of the ridge distribution at low pt- 



V. CONCLUSION 

We have given an interpretation of the ridge phe- 
nomenon in pp collisions in terms of soft partons on which 
very little is known. By drawing on what we do know 
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about the soft partons in nuclear collisions, we are led 
to the implication that a dense medium can be created 
even in pp collisions at 7 TeV. The primary input in our 
approach to explaining the observed ridge yield is the 
assertion that the correlation is of the factorizable form 
R(l)R(2), where R(i) is the response of the zth soft parti- 
cle to the unobserved jet, so that two independent trans- 
verse correlations of the semihard-soft type can lead to a 
net soft-soft correlation in C?2(l,2). 

The success of our approach applied to pp collisions at 
7 TeV suggests that (a) the medium can be responsive to 
semihard jets, (b) there can be azimuthal anisotropy, (c) 
the pr spectrum in the ridge is harder than that of the 
inclusive, and (d) that hadronization is by recombination. 
None of the above rely on the validity of hydrodynamics 
for pp collisions, or the existence of intrinsic long-range 
longitudinal correlation, and all of them can be checked 
by further experimental measurements. The last item 
cannot be checked directly, but one of its consequences 
is that the p/ir ratio can be large, which is a property 
of all recombination/coalescence models [3{|. We expect 
the p/ir ratio in the ridge to increase with px at low 
Pt in pp collisions at 7 TeV, although the rate of that 
increase depends on the soft parton density, on which 



we have insufficient knowledge to predict. A ratio larger 
than 0.2 cannot be explained by fragmentation. Thus 
the experimental determination of the p/it ratio in the 
ridge will be very interesting and should provide further 
insight on the structure and origin of the ridge. 

The basic issue that the observation of a ridge by CMS 
has opened up is whether a system of high density soft 
partons can be created in pp collisions. The system 
may be too small for the applicability of hydrodynam- 
ics, but azimuthal anisotropy can nevertheless exist for 
small systems in non-central collisions, so consequences 
on (f> asymmetry should be measurable, as the ridge struc- 
ture on the near side demonstrates. Our consideration of 
ridge formation as being generated by semihard jets ap- 
plies to both hadronic and nuclear collisions. Thus we 
go further to suggest that even in single-particle distri- 
bution in pp collisions at LHC there may exist a ridge 
component that contains all the <p dependence, as found 
in Au-Au collisions [U H, HI- 

This work was supported in part, by the U. S. Depart- 
ment of Energy under Grant No. DE-FG02-96ER40972 
and by the National Natural Science Foundation of China 
under Grant No. 10775057 and 11075061. 



CMS Collaboration, J. High Energy Physics 09, 091 
(2010). [20 
J. Adams et al, (STAR Collaboration), Phys. Rev. C 73, 
064907 (2006). [21 
B. I. Abelev et al., (STAR Collaboration), Phys. Rev. C 
80, 064912 (2009). [22 
B. Aver et al, (PHOBOS Collaboration), Phys. Rev. [23 
Lett. 104, 062301 (2010). [24 
B. I. Abelev et al, (STAR Collaboration), Phys. Rev. 
Lett. 105, 022301 (2010). [25 
M. Daugherity (for STAR Collaboration), J. Phys. G: 
Nucl. Part. Phys. 35, 104090 (2008). [26 
E. Shurvak, larXiv: 1009.46351 
A. Dumit ru et aZ.,[ar Xiv: 1009.5295 [27" 
P. Bozek, l arXiv: 1010.04051 

J. Putschke, J. Phys. G: Nucl. Part. Phys. 34, S679 [28 
(2007). 

R. C. Hwa, in Quark-Gluon Plasma 4, edited by R. C. [29 
Hwa and X. N. Wang (World Scientific, Singapore, 2010). [30 
S. Gavin, L. McLerran, and G. Moschelli, Phys. Rev. C [31 
79, 051902(R) (2009). 

A. Dumitru, F. Gelis, L. McLerran, and R. Venugopalan, [32 
Nucl. Phys. A 810, 91 (2008). [33 
K. Dusling, F. Gelis, T. Lappi and R. Venugopalan, Nucl. 
Phys. A 836, 159 (2010). [34 

B. Alver and G. Roland, Phys. Rev. C 81, 054905 (2010). 
R. C. Hwa and C. B. Yang, Phys. Rev. C 70, 024905 [35 
(2004). 

C. B. Chiu, R. C. Hwa and C. B. Yang, Phys. Rev. C 78, [36 
044903 (2008). 

G. Aad et al, (ATLAS Collaboration), Phys. Lett. B [37 
688, 21 (2010). [38 
V. Khachatryan et al., (CMS Collaboration), Phys. Rev. [39 



(STAR Collaboration), 
(STAR Collaboration), 



Lett. 105, 022002 (2010). 

K. Aamodt et al, (ALICE Collaboration), Phys. Lett. B 
693, 53 (2010). 

C. -Y. Wong, Phys. Rev. C 78, 064905 (2008); 80, 034908 
(2009). 

G. Moschelli and S. Gavin, Nucl. Phys. A 836, 43 (2010). 
K. Werner et al, Phys. Rev. C 82, 044904 (2010). 

A. Feng (for STAR Collaboration), J. Phys. G: Nucl. 
Part. Phys. 35, 104082 (2008). 

D. Netrakanti (for STAR Collaboration), Nucl. Phys. A 
830, 681c (2009). 
M. M. Aggarwal et al.. 
larXiv: 1004.23771 

H. Agakishiev et al, 
larXiv:1010.0690l 

C. B. Chiu and R. C. Hwa, Phys. Rev. C 79, 034901 
(2009). 

R. C. Hwa and L. Zhu, Phys. Rev. C 81 , 034904 (2010). 
R. C. Hwa and L. Zhu. larXiv: 1101. 13341 

B. B. Back et al. (PHOBOS Collaboration), Phys. Rev. 
Lett. 91, 052303 (2003). 

C. B. Chiu and R. C. Hwa. larXiv: 1012.3486 1 

J. Adams et al, (STAR Collaboration), Phys. Rev. Lett. 
95, 152301 (2005). 

C. B. Chiu and R. C. Hwa, Phys. Rev. C 72, 034903 
(2005). 

R. P. Feynman, Photon- Hadron Interaction (Benjamin 
Press, Reading, MA, 1972). 

J. D. Bjorken, SLAC Summer School, SLAC-167 vol. I 
(1973) (unpublished). 

R. C. Hwa, Phys. Lett. B 666, 228 (2008). 

K. P. Das and R. C. Hwa, Phys. Lett. 68B, 459 (1977). 

R. C. Hwa and C. B. Yang, Phys. Rev. C 67, 034902 



6 



(2003); V. Greco, C. M. and P. Levai, Phys. Rev. Lett. S. A. Bass, Phys. Rev. Lett. 90, 202303 (2003). 

90, 202302 (2003); R. J. Fries, B. Miiller, C. Nonaka, and 



